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Abstract 

A framework combining Yang-Mills dynamics of the pre-equilibrium glasma with relativistic 
viscous hydrodynamic evolution of the quark-gluon plasma and hadron gas phases is presented. 
Event-by-event fluctuations of nucleon positions and color charges are taken into account, lead- 
ing to negative binomial fluctuations of gluon multiplicities. Experimental anisotropic flow coef- 
ficients V2-V5 of charged hadron distributions in heavy-ion collisions at the Large Hadron Collider 
i-Sh 1 are well described. Furthermore, event-by-event distributions of Vz, V3 and V4 measured by the 

ATLAS collaboration are reproduced. 
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1. Introduction 

Heavy-ion collisions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron 
Collider (LHC) allow for systematic exploration of the high temperature many-body dynamics 
of a non-Abelian quantum field theory. The study of multiplicity fluctuations and anisotropic flow 
harmonics v„ at both RHIC HI 01 and LHC yl 51 IH potentially allows to obtain information on 
the strongly correlated non-equilibrium glasma regime and the transport properties of the nearly 
equilibrated quark-gluon plasma and hadron gas phases. To do so, a sophisticated theoretical 
description of all stages of the evolution of the complex system created in heavy-ion collisions is 
needed. We present a promising theoretical framework for this task with the combination of the 
IP-Glasma model with relativistic viscous hydrodynamics. 

2. IP-Glasma model 

The IP-Glasma model combines the IP-Sat (Impact Parameter Saturation) model 
of high energy nucleon and nuclear wavefunctions with the classical Y ang -Mills (CYM) dynam- 
ics of the glasma fields produced after the heavy-ion collision US [H [H 51. It relates the 



deeply inelastic scattering (DIS) constrained nuclear dipole cross-sections to the initial classical 
dynamics of highly occupied gluon "glasma" fields after the nuclear collision. Given an initial 
distribution of color charges in the high energy nuclear wavefunctions, the IP-Glasma framework 
computes the strong early time multiple scattering of gluon fields by event-by-event solutions of 
Yang-Mills equations. 

The IP-Glasma model naturally includes the effect of several sources of quantum fluctua- 
tions that can influence hydrodynamic flow: fluctuating distributions of nucleons in the nuclear 
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Figure I: Probability distributions of gluon multiplicities dN g /dy. Left: RHIC. Shown are also the distributions for some 
limited ranges of N pMt , which are described by negative binomial distributions. Right: LHC with centrality classes. 



wavefunctions and intrinsic fluctuations of the color charge distribution. This results in "lumpy" 
transverse projections of the gluon field configurations that vary event to event. The scale of this 
lumpiness is given on average by the inverse nuclear saturation scale Q s which corresponds to 
distance scales smaller than the nucleon size. The details of the implementation can be found in 

3. Multiplicities 

In Fig. Q] we present the probability distribution of dN g /dy at RHIC and LHC energies. An 
essential ingredient is the probability distribution of impact parameters, which is determined by 
the Glauber model. One could in principle compute this distribution in the Glasma framework, 
but a first principles computation is extremely difficult (see J3l). For RHIC we also show five 
distributions obtained by constraining the number of participants to demonstrate that their shape 
resembles negative binomial distributions. For LHC we indicate ranges for different centrality 
classes. 

4. Transition to hydrodynamics 

When we switch from the CYM description to hydrodynamics we construct the fluid's initial 
energy momentum tensor T^-. — (e + P)u fi u v - Vg^ + TP" from the energy density in the fluid's 
rest frame s, the flow velocity u 11 , and, using an equation of state, the local pressure f at each 
transverse position, s and are obtained by solving M/jT'cym = suV, usm § tne f act tnat u>1 K a 
time-like eigenvector of T^ YW and satisfies u 2 — 1 . For the present study we do not extract the 
very non-equilibrium IW but assume an efficient instability driven isotropization of the system 
and set the initial YV lv to zero. A 3+1 dimensional Yang -Mills simulation including quantum 
fluctuations can hopefully provide the full dynamics for this mechanism in the future. Further 
details on the hydrodynamic simulation can be found in II 1511 . 

We determine centrality classes using the multiplicity distribution of gluons (Fig. Q] right) 
much alike the procedure followed by the heavy-ion experiments using charged particle multi- 
plicity distributions. 
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The hydrodynamic sta ge, including a Cooper-Frye freeze-out and resonance decays, is sim- 
ulated using music Using the final particle distributions we determine 



v k (Pt) = <cos(n(0 - ^„))) (1) 
in every event using the exact event plane 

1 <sin(«0)> 

if/„ = - arctan — — . (2) 

n (cos(n0)) 

where (■) is the average over the smooth particle distribution function. 



5. Flow results 

Calculated root-mean-square (rms) V„(pt) for 10-20% central and more peripheral collisions 
agree very well B15I1 with experimental data from the ATLAS collaboration [4], which are deter- 



mined using the experimental event-plane method that agrees well with the rms values 1 20, 21]. 
Integrated rms v>2, V3 and v\ also agree very well 11511 with available v„{2] results (obtained using 
two-particle correlations) from the ALICE collaboration yfl. 

Here, we present event-by-event distributions of V2, V3 and V4 for two centrality classes com- 
pared to results from the ATLAS collaboration 1 2(J 21 ] in Fig. [2] Distributions are scaled by their 



mean to be able to compare to initial eccentricity distributions at the same time. We find that the 
initial eccentricity distributions are a good approximation to the distribution of experimental v n . 
Exact agreement with the data in the 20-25% central bin, however, is only achieved with the v„ 
distributions obtained after hydrodynamic evolution. 



6. Summary 

The IP-Glasma+Music model produces anisotropic flow in excellent agreement with experi- 
mental data from the LHC. It is particularly remarkable that the event-by-event distributions of 
v„ agree with experimental results from the ATLAS collaboration, indicating that the fluctuations 
in the initial state are described accurately by the model. The successful description of a wide 
range of data within this model provides a framework to precisely determine key aspects of the 
complex dynamics of heavy ion collisions. 
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Figure 2: Scaled distributions of vi, V3 and v 4 as well as £2, £3 and e 4 compared to experimental data from the ATLAS 
collaboration ['20, 21]. Using 750 (0-5%) and 1300 (20-25%) events. Bands are systematic experimental errors. 
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